function para_gms = gen_gms_perturb(para_pa_init, para_rand, N_perturb, objga)

[Tn, cv] = objga(para_pa_init);
para_k = para_pa_init;
if Tn<cv
    para_gms = para_k;
else
    para_gms = [];
end

for k = 1 : N_perturb
    para_k2 = para_k + para_rand(k, :);
    [Tn2, cv2] = objga(para_k2);

    if Tn2<cv2
        para_k = para_k2;        
        para_gms = [para_gms; para_k2];
    end
end